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Abstract 

We consider a simple model of large- N c QCD defined by a spectrum consisting of an infinite set 
of equally spaced zero-width vector resonances. This model is an excellent theoretical laboratory 
for investigating certain approximation schemes which have been used recently in calculations of 
hadronic parameters, such as the Minimal Hadronic Approximation. We also comment on some 
of the questions concerning issues of local duality versus global duality and finite-energy sum rules. 


1 Introduction 


The suggestion to keep the number of colours N c in QCD as a free parameter was made by’t Hooft jlj 
as a possible way to approach the study of non- perturbative phenomena. In spite of the efforts of 
many illustrious theorists who have worked on the subject, QCD in the large-iV c limit still remains 
unsolved. However, many interesting properties have been proved^ which suggest that, indeed, the 
theory in this limit has the bulk of the non-perturbative properties of three-colour QCD. In particular, 
it has been shown that, if confinement persists in this limit, there is spontaneous chiral symmetry 
breaking Q. 

The spectrum of the theory in the large-iV c limit consists of an infinite number of narrow stable 
meson states which are flavour nonets This spectrum, however, looks a priori not quite like the one 
in the real world. The vector and axial-vector spectral functions measured in e + e~ —» hadrons and in 
hadronic r decay show indeed a richer structure than just a sum of narrow states. There are, however, 
many instances where one is only interested in observables which are given by weighted integrals of 
hadronic spectral functions. In these cases, it may be enough to know a few global properties of the 
hadronic spectrum in order to have a good interpolation. Typical examples of that are the coupling 
constants of the effective chiral Lagrangian of QCD at low energies, as well as the coupling constants 
of the effective chiral Lagrangian of the electroweak interactions of pseudo-scalar particles in the 
Standard Model. It has been shown in some examples i 0-! 10 that inserting the hadronic 
spectrum of large-iV c QCD as an approximation to the real hadronic spectrum provides rather good 
results. It is in this sense that large-lV c QCD seems to be a very useful phenomenological approach 
for understanding non-perturbative QCD physics at low energies. 

In most cases of interest, the Green’s functions which govern the low-energy constants of the chiral 
Lagrangian are two-point functions with zero-momentum insertions of vector, axial-vector, scalar 
and pseudo-scalar currents. The higher the power in the chiral expansion, the higher the number 
of insertions. Furthermore, they are order parameters of spontaneous chiral symmetry breaking; i.e. 
they vanish, in the chiral limit where the light quark masses are set to zero, order by order in the 
perturbative vacuum of QCD. That implies that they have a power-like fall-off in l/Q 2 at large Q 2 ; 
where Q 2 = — q 2 > 0 in the euclidean region, with q the four-momentum flowing through the Green’s 
function. When N c —> oo, within a finite radius in the complex Q 2 plane centered at the origin, 
these Green’s functions only have a finite number of poles. The restriction to the minimal number 
of poles required to satisfy the known short-distance and long-distance QCD constraints results in a 
simple approximation which we call the minimal hadronic approximation (MHA) to large- N c QCD. 
In the cases which we have been able to test, this minimal approximation gives already very good 
phenomenological results. The MHA is in principle improvable: considering more terms in the OPE 
of the two currents in the underlying Green’s function through which the g-momentum flows provides 
extra algebraic sum rules which can be used to fix the extra hadronic parameters. Also, knowing more 
terms at small-Q 2 values, in the chiral expansion of the same underlying Green’s function, provides 
extra algebraic sum rules of a similar nature. Unfortunately, in most cases of interest, it is difficult to 
carry on these improvements in practice, and the nature of the MHA in real QCD remains unclear. 
This is the reason to resort to simple models, where the full large-iV c -like spectrum is fixed a priori 
and one can study approximations and their improvements explicitly. In this paper we shall study a 
model defined by a spectrum of an infinity of zero-width vector particles whose masses (squared) are 
equally spaced. We consider this spectrum quite realistic. As a matter of fact, for large masses, this 
is the spectrum of QCD at large N c in 2 dimensions (where the theory can be solved |ll|) and, in 
the case of our four-dimensional world, this is also the spectrum of the daughter trajectories of Regge 
Theory; which has long been suspected to be a good effective description of QCD in the large-W; 
limit. 

To the best of our knowledge, this type of model was first proposed in the context of QCD in 
ref. |l2|]. More recent studies confronting phenomenology can be found in ref. |l^, [L4|. Also recently 
Shifman p5[ has emphasized the usefulness of this model for investigating the question of local duality 
versus global duality. The model has very intriguing connections with functions which appear profusely 
in analytic number theory. 

1 For recent reviews, see e.g. the book in ref. O] and the lectures in ref. [H. 
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The paper is organized as follows. In section 2, we introduce the model with equally spaced, 
infinitely narrow resonances. We then take the vector-current two-point function (more precisely, the 
Adler function) as a laboratory, and calculate it in this model, including its OPE and chiral expansions. 
In section 3 we consider a specific form of the minimal hadronic approximation. After some discussion 
of its nature, we give explicit numerical examples of how the MHA approximates the Adler function 
of our model. Finally, in section 4, we compare the hadronic branching ratio for r decay between the 
model, the OPE and the MHA, in order to see how the MHA performs with respect to local duality 
and finite-energy sum rules. 


2 The model 

Let us consider the vector-vector correlator defined by 

* J d 4 xe iq ' x (0\T r(0)) |0> - (gV - g^q 2 )H(Q 2 = -q 2 ), (2.1) 

with the electromagnetic current of light quarks, for example. The model for this correlator consists 
of an infinite number of equally spaced narrow states as illustrated in Fig. 1 below. With Mq the 
mass of the lowest state and <r the spacing mass, the spectrum is given by 

1 

—Imll(t) — A a 2 S(t — Mq — na 2 ) , (2.2) 

7r z ' 

n—0 

with A an arbitrary normalization factor that will be determined below. 



Fig. 1 The vector spectral function in the toy model of large-N c QCD. 


Of special interest for our discussion is the Adler function defined as follows: 
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where £ ^2, — ^ is the so-called generalized Riemann zeta function sometimes called the 

Hurwitz functior^. We can choose the normalization factor A so that the large-Q 2 behaviour of the 
Adler function reproduces the parton like asymptotic freedom behaviour predicted by perturbative 
QCD. This can be readily seen using the Euler- Maclaurin formula: 


1 1 1 



and the integral result 
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MS 
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+ N 



(2.5) 


( 2 . 6 ) 


(2.7) 


Letting TV —> oo we get the result 

lim A(Q 2 )^A, (2.8) 

Q 2 —> oo 

which, for A = reproduces the asymptotic parton limit result. For simplicity, we shall take 

A = 1 from here onwards. 

There is an integral representation of the Hurwitz function which will be very useful for our 
purposes: 

1 t S ~ 1 

C(s, z) = . / dte~ zt - -T (Re s > 1 and Re z > 0). (2.9) 

T(s) J 0 1 - e 


2.1 Operator Product Expansion of the model 


In QCD, the behaviour of the Adler function in the physical vacuum, for large euclidean Q 2 values, 
can be obtained from the operator product expansion (OPE) of the two currents in Eq. (EH) and it 
results in an expansion in 1/Q 2 powers a la Shifman, Vainshtein, Zakharov fl7|. In the model this 


expansion follows readily from the integral representation (2.6) applied to ( ^2 


q 2 +m 2 


Using the fact that 


±^ = J2B n (x)-z\ 
— 1 n\ 


n —0 


is the generating function of the Bernoulli polynomials ,[] we have 


o b 


= e 




Ar-EM 1 - 


n—0 


M 2 \ 1 


- t n 


( 2 . 10 ) 


( 2 . 11 ) 


and using the symmetry property in Eq. (S plus the Laplace transform integral representation 

1 1 


s n+1 r(n + l)7 0 


( 2 . 12 ) 


we get the required OP E-like expansion 

2 It also coincides with 'ip '), where 'ip'(z) = logr( 2 :) and T( 2 :) is Euler’s Gamma function. 

3 For the sake of completeness we list a few useful properties of the Bernoulli polynomials in an Appendix. 
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(2.13) 
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The explicit expansion for the first few terms is then 
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In QCD, and in the chiral limit where the light quark masses are set to zero, there is no l/Q 2 term 
in the OPE, because there is no local operator of dimension d = 2. We can then make the model look 
like QCD more closely by imposing this property in Eq. ([2.14) as well, which requires that 


Mn = -< 


(2.15) 


The model has then only one mass scale a 2 ; furthermore, because of the property of the Bernoulli 
polynomial that 

B 2 „+i( 1/2) = 0, for n > 0, (2.16) 

the equivalent of the OPE of the Adler function in the model, has only even 1/Q 2 powers: 


l/cr 2 \ 2 7 /er 2 \ 4 31 (a 


12 \Q 2 


A (oo) (Q 2 ) 1 - ttt ( Tjo ) + E -7^7 E 


240 \Q 2 


1344 \Q 2 


2 \ 6 


(2.17) 


and the vacuum expectation values appearing in the OPE are then given by the values of the Bernoulli 
polynomials at x = 1/2. For the rest of the paper we shall keep Mg = cr 2 /2 and comment on the 

relevance of this constraint whenever necessary. _ 

It is easy to see that the OPE series of the model in question is a divergent series. From Eq. ( A.15 ) 
it follows that 


£Ml/2) = (—l) iV+i 2(2iV)! 
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(2tt) 2N *—' n 

' ' n =1 
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(-i) r 


2N 


and the sum of the series on the right-hand side is known: 


^ (-1)" (l _ 2 2Ar - 1 ) 7 r 2Ar 
zL ~^N~ - -77777-™ | , 


(21V)! 


which results in 


B 2N (l/2) = {-\) N [\-^-^j\B 2N \, N >0. 


From this, it is easy to see that the series is divergent because 

2JV 
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where we have used the property that asymptotically, 
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At orders 
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the OPE ceases to be an improving expansion. 

The OPE series is divergent, but summable. In fact, the Borel sum of the OPE series is precisely 
the Adler function defined by the generalized Riemann zeta function. To see this, recall that the Borel 
sum of the series in Eq. ( 2.13|) is defined as follows 


£(-l)"B n (l/2) 

n —0 


a 

Q 2 


dte 'E 1-1 -) 1 


n—0 


i B n (l/2) (a 2 _ 

Q 2 


(2.24) 


If now, we use again the symmetry property in Eq. (A.13) and the second of Eq. (2.11); plus the 

~ 2 - - 

change of variables t = one obtains the desired result. 


2.2 Chiral Perturbation Theory expansion of the model 


We can also find, in this model, the equivalent of the chiral perturbation theory (yPT) expansion of 
the Adler function; i.e., the Taylor expansion at small Q 2 values. For that, it is sufficient to expand 

_ q 2 t —-n 

the term e ^ in the integral representation of Eq. (|2.9|) with the result 


71—0 


- 4 ( 0 ) (Q 2 ) = E(— !)" (%) (n+l) tln + 2, 1 -^-) . 


2 \ n + 1 


Ml 


With Mq = icr 2 , and using the property that 


C( s , 1/2) = (2 s — 1) C( s ) > 


for 


s > 1, 


the series becomes 


71=0 


A {0 \Q 2 ) = E(- 1 )" ( n + 1 ) 2 (n+2) -:L C(n + 2)[% 


2 \ ™+l 


(2.25) 


(2.26) 


(2.27) 


i.e., the low-energy constants of the Adler function are proportional to Riemann’s zeta functions of 
increasing argument. The explicit chiral expansion for the first few terms gives 


■4 <0> (Q 2 )-4(^)-14<(3)(^t -124«5)(^) + G [^ 
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Numerically, 


A {0 \Q 2 ) -> 4.9348 (0^ - 16.8288 ^ 


48.7045 


Q‘ 


- 128.579 


Q 5 


(2.28) 


(2.29) 


q2 

The chiral expansion is convergent within a radius ^ 

It is indeed amusing that a model of a well defined limit of a rather complicated quantum field 
theory like QCD, has such direct connections with the basic functions of analytic number theory and 
transcendental numbers. 


3 The Minimal Hadronic Approximation. 


We want to study how good is the approximation which consists in replacing the exact spectral 
function in Eq. (2.2) by a finite set of low states, say m+ 1, with arbitrary masses A/ 2 , and lumping 
the rest into the equivalent of a perturbative QCD continuum starting at a mass squared value Sq: 


—ImII(f) 

7 T 


(m+ 1) 

MHA 


= a 2 E <*(* - M 2 ) + 0(t - s 0 ) • 

n =0 


(3.30) 
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The approximated spectrum is illustrated in Fig. 2. We call this the minimal hadronic approximation 
(MHA), by analogy to the approximations which we have been making in real large- N c QCDjf] 



Fig. 2 The Minimal Hadronic Approximation (MHA) to the vector spectral function in the toy 
model of large-N c QCD. The successive approximations correspond to one state plus continuum (the 
first blue band), two states plus continuum (the second blue band), three states plus continuum, etc. 


The MHA spectrum leads to the following form for the Adler function: 


A(Q 2 ) 


(ra+1) 

MHA 


Q 2 y 1 

he (K + ML 


Q 2 

Q 2 + So 


(3.31) 


where the second term is the one induced by the perturbative continuum. The OPE-like expansion 
of the MHA approximation to the toy model can be easily obtained from the expression above, with 
the result 


A(~>(Q 2 )| mha = !+£(-!)" 

N=1 



m 


n —0 




(3.32) 


We fix the onset of the continuum so to match the leading l/Q 2 term of the OPE for the Adler 
function, which vanishes. This results in the condition^ 


s 0 = (m + l)tr 2 


(3.33) 


The masses Mfj of the explicit narrow states are then fixed so as to satisfy more and more 
constraints provided by the OPE and the \PT expansions of the underlying theory. This results in a 
system of equations: 


{m + l) N 


m 


n =0 



N-l 

= B n { 1/2), 


for N = 2,3,..., 


(3.34) 


4 See ref. ]l8| for a r ecent review, and refs. 0 to Jic| ] for details. 

5 If Mq / 1/2, Eq. ( |3.33| ) reads so = cr 2 (m + 1/2 + M 0 2 ). Notice that so should always be larger than the mass of 
the last resonance kept explicitly in the spectrum. 
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from the OPE constraints; and another system of equations 


N(m + 1) N 


, ( a ^_\ N+1 
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+ ^ V M 2 ) 
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(3.35) 


n —0 


- for N = 1 

(n + !)*+!’ 


from the %PT constraints. The idea is to solve these equations for successive values of the number 


m of explicit narrow states. Looking at Eq. (3.35) one sees that the first term on the left- hand side 
vanishes in the m —► oo limit. It is clear that, in this limit, the system of equations ( |3.35 ) has as a 
solution M 2 = cr 2 (n + 1/2), i.e. the true spectrum. 

We also find (numerically, including up to nine resonances) that wh en on e substitutes the masses 
found from only the yPT constraints into the left-hand side of eq. ( |3.34 ), the resulting numbers 
converge to the values given by the right-hand side (i.e. the B^{l/2)). This implies that in principle, 
the MHA will converge to the exact result for all values of Q 2 , even if one uses only small Q 2 
information. However, the convergence is rather slow for large Q 2 . Therefore, in practice, one would 
expect better convergence if one also uses constraints coming from the OPE in determining the 
parameters M 2 of the MHA. 

It is less clear what the nature is of the system of equations following from using only OPE 
constraints since, unlike the \PT expansion, the O PE is an asymptotic expansion (see section (2.1)). 
For a given number of OPE constraints, 2 N, Eq. ( 2.23| ) gives a value of Q 2 below which the MHA 
obtained by solving the set of equations ( [3.34 ) may be wrong J^j Conversely, gi ven a value of Q 2 above 
which one wishes to obtain a good approximation to the exact model, eq. (2.23) gives an estimate 
of the maximum useful order of the OPE. However, as already alluded to above, we find it quite 
plausible to expect that using a few constraints coming from the OPE in addition to those from yPT 
should speed up the convergence in the region of large Q 2 , relative to a MHA solely based on yPT 
constraints. In practice, quantities which typically one encounters in QCD correspond to integrals of 
the Adler function weighted by a known function, and inspection of the weight function should dictate 
the optimal choice of short-distance constraints versus long-distance constraints. For illustration 
purposes, in the examples which follow below, we shall always combine both short- and long-distance 
constraints. An important comment, however, is that when one uses only chiral constraints, the masses 
M n appear to converge monotonously to the exact values, while if one also uses OPE constraints, this 
appears not to be true. Even so, the approximate Adler function, in all the numerical experiments we 
have examined, always gets better for all values of Q 2 if one includes more resonances in a systematic 
way. 

We shall now construct explicitly the leading and next-to-leading Minimal Hadronic Approximation 
to the model. 


3.1 One State plus Continuum 
The onset of the continuum is now 

and the first OPE constraint becomes 


The numerical value of the solution is 


s 0 = cr 


M 2 

1-2-f 

a z 


1 

12 


Mq = 0.541667er 2 , 


(3.36) 


(3.37) 


(3.38) 

6 Although not necessarily so. The Stieltjes function J^° dt -^—, which has 0 ( — l) n n\z n as the asymptotic 

series for 2 : —*■ 0, can actually be approximated in a convergent sequence of steps with this method |19l. 
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not quite the exact value Mg = ^a 2 of the real spectrum but not bad at all! In fact, as shown in 
Fig. 3, the simple MHA of one state plus continuum, already reproduces rather well the shape of 
the exact Adler function. For comparison, if M§ is fixed using the first X'PT constraint, one finds 
Mq = 0.504125er 2 . While this is much closer to the exact value Mg = 0.5a 2 , the approximate Adler 
function turns out to be a worse approximation for Q 2 > 1.8a 2 . 



Fig. 3 The exact Adler function of the toy model (the continuous red curve), and the Adler 
function predicted by the MHA with one state plus continuum (the dashed blue curve), are plotted 
versus Q 2 /a 2 . 


3.2 Two States plus Continuum 

We can now tune the masses Mq and Mi to reproduce the leading and next-to-leading Ag terms 
of the underlying OPE. The constraining equations are now 


4-2 


J a 2 " a 2 

n 0 M n 4 n Mf 

^ 8 + 3-IL+3-A 


a’ 


-1 
~L2 ’ 


(3.39) 

(3.40) 


They lead to the solutions: 


M 0 2 = 0.481174a 2 and M 2 = 1.56049a 2 , (3.41) 

instead of the physical values Mg = ^a 2 and M 2 = |a 2 . Notice that the mass of the ground state in 
the two state approximation gets closer to the physical mass. The corresponding Adler function also 
approaches better the exact Adler function. In order to quantify this improvement we plot in Fig. 4 




the difference of functions 


(i) 


( 2 ) 


A« = A{Q 2 ) - A(Q 2 ) and A (2 ) = A(Q 2 ) - A(Q 2 ) 


versus Q 2 


(3.42) 



0 0.2 0.4 0.6 0.8 1 1.2 1.4 


Fig. 4 The difference of Adler functions defined in Eq. ( | S.fA ) corresponding to the MHA with one 
state plus continuum (the dashed green curve) and with two states plus continuum (the continuous 
blue curve) are plotted versus Q 2 /a 2 . Notice the vertical scale in the plot. 


Another alternative is to use as constraints the leading Aj behaviour of the OPE and the leading 
Q 2 behaviour of the chiral expansion. The constraining equations are then 


Mg M? 
4 - 2—- - 2—- 


a“ 


<j 


A 4 

a a 

~M§ + ~Mf 


-1 
12 ’ 
7r 2 

T ’ 


and the corresponding numerical solutions for the masses are then: 

and M 2 = 1.54257cr 2 . 


M 0 2 = 0.499093ct 2 


(3.43) 

(3.44) 

(3.45) 


They approach, even better, the exact masses than those previously obtained in Eq. ( 3.41 ). In fact, 
combining the leading short-distance constraint with the leading long-distance constraint turns out 
to produce an approximated Adler function much closer to the exact one. We illustrate this in Fig. 5 
where the dashed blue curve is the same as the continuous blue curve in Fig. 4 (though plotted on a 
different scale), while the continuous red curve is the one obtained solving the constraints of Eq. ( 3.43 ). 
The improvement is rather remarkable. We would like to emphasize that in the applications of the 
MHA to large- N c QCD that we have made so far in refs. Q to 0 , we have used equivalent constraints 
to the ones discussed in this example. 
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In Fig. 5 we use the following notation: Aindicates the difference between the exact Adler 
function and the one obtained using the MHA with m states plus continuum, with the masses obtained 
by solving a system of equations which consists of i short-distance constraints and j long distance 
constraints; i.e., 


A^ = A(Q 2 ) - A(Q 2 ) 


(m) 


with i SD-constraints 
and j LD-constraints 


i+j = m. 


(3.46) 



Fig. 5 The difference of Adler functions corresponding to the MHA with two states plus continuum. 
The dashed blue curve is the same as the one in Fig. f, though on a different scale, while the continuous 
red curve is the one obtained using as input masses the solutions in Eq. (3-45). 


We have constructed the MHA with up to 11 resonances. With a similar number of short-distance 
and long-distance constraints, the maximum difference between MHA and the exact Adler function is 
always decreasing by half an order of magnitude (in average) every time a new resonance is explicitly 
included in the MHA. This maximum difference occurs for Q 2 at roughly 1/3 of the number of 
resonances included. For instance for 11 resonances this maximum difference is of the order of 10~ 8 
and happens at Q 2 /a 2 — 4. 


4 Finite—Energy Sum Rules and Hadronic Duality. 

In every day life QCD, the paradigm of a finite-energy sum rule (FESR) is the one provided by the 
hadronic decay of the t, which relates the branching ratio 

K, = r(r- ^ -r Hadrons) 
i(t —> v T e v e ) 
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to the vector and axial-vector hadronic spectral functions. In the chiral limit, this relation is simply 
given by the integral]] 


dt 


K T (s)= — l-~ 1 + 2- —Imn(t), 


t\ 1, 


(4.48) 

for s = to 2 . Sum rules of this type have been extensively studied in the literature^ and proved 
very useful for extracting the value of the QCD running coupling constant. The toy model offers an 
interesting theoretical laboratory to check these sum rules. Some observations in this respect have 
already been made by Shifman [p~5| . For the sake of clarity, we reproduce part of his analysis in our 
discussion. 


Equation (4.48) is a combination of moment sum rules of the type 

M M (s)= [ -f-V-Imn(t). 

Jo s W k 


(4.49) 


In the toy model, with the input spectral function in Eq. (2.2), these moments can be explicitly 
calculated with the result 


M^ p \s) = a 2 


(n + i) p a 2 P 


CP + 1 


n =0 


where 0 < x < 1 is the fractional part of -%■ — i; in other words, 


s = (m + x + -)a 2 , for m = 0,1,2,- 


The exact lZ T (s) in the toy model of large-7V c is then 


(4.50) 


(4.51) 


’W’-HOr ' 


1 f fj 2 x 4 

x(l - x)(l - 2x) + - l — 


x 2 (l - x) 2 - 

v ’ 16 


(4.52) 


This expression could be viewed as the exact lZ T (s) calculated from the data in an imaginary world 
in which the r mass may be dialed to be as large as one wishes. This exact 7 Z T (s) is plotted in 
Fig. 6 versus the ratio s/c r 2 (the red continuous curve ) and as it can be seen, it sh ows a damped 
oscillatory behaviour induced by the convolution in the integrand of 7 Z T (s) in Eq. ( 4.48| ) of a decreasing 
polynomial in s and an infinite set of equally-spaced “kicks” (i.e. the narrow resonances). Should the 
number of resonances stop at a certain finite number to, beyond which the continuum takes over, the 
oscillations for > to + ^ would disappear and turn into an ordinary power fall-off towards the value 
of the parton model, i.e. 1/2. 

It is interesting to compare the exact 7Z T (s) function in Eq. ( 4.52| ) with the one predicted by the 
OPE of the same model. Emulating what is being done in the real-world analysis of r decay, this 
comparison can be obtained from an integral representation in the complex z-plane (where z = Q 2 /s), 


T^t{s) = —— ® dz (1 + z) 2 (1 — 2 z) n(zs) 

J\z 1=1 


1 

27ri 


'1*1= 

dz 

1 * 1=1 z 


- z ■ 

2 


2, Z 


A(zs) 


(4.53) 

(4.54) 


by inserting the OPE for the Adler function in Eq. m- Doing this one notices that, besides the 
partonic contribution, only the powers proportional to 1 /Q 2 1 1 /Q 6 and 1 /Q * 8 in the OPE of the Adler 
function contribute to the integral, with the result 

' Here we only consider the vector, and not the axial—vector, part. 

8 See e.g. refs. pi| pit. 
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(4.55) 


n 


OPE 

T 



7 

480 



The curve corresponding to this TZ GPE (s) function is also plotted in Fig. 6 (the blue dashed curve). 
It interpolates rather well the exact curve but it fails to reproduce its oscillations Jl5| . This is another 
manifestation of the violation of local duality (i.e. point-wise; as opposed to global duality, i.e. on 
average) between the quark -gluon picture represented by the OPE and the hadron picture represented 
by the full spectrum- based solution. 



Fig. 6 Plot of the branching ratio lZ T (s) defined in Eq. ( tf.JSl ). The continuous red curve is the 
one predicted by the exact toy model; the dashed blue curve is the one predicted using the OPE. 


However, one should also notice in this figure that there are “optimal” points at s/a 2 ss m/ 2, m = 
2,3,4,... where the difference between the oscillation and the OPE curve nearly vanishes. These 
“optimal” points are very useful because they allow the OPE to produce a very accurate result even 
at very low values of s/a 2 ! We shall now elaborate on these special points and explain where they 
come from. To this end it is useful to identify clearly from 1Z T the combination which goes into 1Z GPE 
and the one which goes into 7?.)? sc . This can be readily obtained by rewriting lZ T (s) in Eq. (4.52) as 
a linear combination of Bernoulli polynomials in the x variable. Using Eqs. (A.8) and (A.9) we can 
rewrite 

Ur(s) = TZ GPE (s) + TZ GSG (s) , (4.56) 


with 7?° pe (s) as given in Eq. ([4.55D, and 


n‘ 


osc 


(«) = - 


t) B3(i) + Ht) Bt(x) 


(4.57) 
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The oscillatory nature of this expression follows from the Fourier decomposition of the Bernoulli 
polynomials, as given in Eqs. (A.15) and (A.16), with the result 


n: 


osc 


m = -i 2 £ T 


sin [27m (£-£)] „ / CT 


(27rn) 


+ 2 — 
s 


cos [27m (4y - 1)] 
(27m) 4 


(4.58) 


where we have used the fact that x, defined in Eq. ( 1.5l| ), is equivalent to 4* — 1 in the argument of 
sin [27rna;] and/or cos [27rnx]. Equation (4.58) shows the oscillations decomposed in terms of harmonics. 
Interestingly, one sees that the leading oscillation at large s, i.e. ~ 1/s 3 , cancels exactly for values of 
s/a 2 which are integers or half-integers, i.e. = -|, A; = 0,1,2, • • • [] The set of half-integers coincide 
with the exact location of the resonances and would correspond in the case of true QCD to the masses 
in the large- N c limit, which is informatio n one does not have in real life. However, the integer values 
coincide with the condition on s in Eq. ( 3.33 ), namely s = (m + 1 )a 2 . This is more useful since it 
only relies on the matching of the MHA onto the first term of the OPE which is always available. Of 
course, it is not known to what extent a similar relation there exists between the leading oscillatory 
term and the l/Q 2 term in the OPE for the case of large- N c QCD; these features may very well be 
generic, and should be kept in mind when dealing with phenomenological studies of finite-energy sum 
rules in QCD. Indeed, the anal ysis of Ref. |^2j shows that the data seem to support this. 

The representation in Eq. (4.58) neatly disentangles the effect of the spectrum (the infinite sum 
over the oscillatory terms) from its projection on the particular observable we are considering; in this 
case 7 Z T , which is only sensitive to the two powers (cr 2 /s) 3 and (er 2 /s) 4 . 

Looking at Fig. 6 makes it clear that an analysis of the data at a given value of s based solely on the 
OPE curve induces an error which is proportional to the size of the nearby oscillation. Consequently 
it is highly desirable to develop an approximation scheme which, at least in principle, can reproduce 
these oscillations. We shall now see that the MHA is such a scheme. 

One example of the prediction from the MHA for 1Z T can be seen in Fig. 7. This example is the 
case of four resonances which are matched onto the first two long-distance and the first two short- 
distance constraints. Since the MHA has been constructed by matching on to the first few terms of 
the OPE, it has to reproduce the corresponding prediction of the OPE for 1Z T as soon as s > so, 
so being the onset of the continuum of the MHA. Since 1Z T can only read the inverse powers up to 
and including l/Q 8 , this means the MHA will reproduce the result of TZ^ FE provided at least three 
resonance masses are matched on to the corresponding three short-distance constraints corresponding 
to the powers l/Q 4 , l/Q 6 , l/Q 8 of the OPE (recall that the l/Q 2 power is always used to determine 
so). Since in the MHA curve in Fig. 7 only the first two inverse powers of Q 2 have been used, this 
explains the small difference between the MHA and the OPE curves in Fig. 7 at s = So = 4<r 2 . One 
can check that this difference disappears when the three short-distance constraints coming from the 
l/Q 4 , l/Q 6 , l/Q 8 powers of the OPE are used. This is not surprising. It is also not surprising that 
the MHA fits so well the low-energy region in Fig. 7, since the MHA has been matched onto the first 
two chiral powers at low energies. What is more interesting is how well the MHA also describes the 
intermediate region of s with all the oscillations. This result clearly goes beyond the input. 

®Notice that for the value x = 1/2, the parton model and the exact result in Eq. 

„ 2 

However, this is no longer true when Mq 7 ^ 1/2. See also footnote 5. 


.52p agree with each other. 
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Fig. 7 Plot of the branching ratio lZ T {s) defined in Eq. ( k^6| ). As in Fig. 5, the continuous red 
curve is the exact result predicted by the toy model, the dashed blue curve is the one predicted using 
the OPE. The black dot-dashed curve is the prediction of the MR A. This curve has been produced with 
four resonances fixed by the first two short-distance and the first two long-distance constraints. 


A useful conclusion, for phenomenological purposes, which emerges from these toy-model exercises 
is that, when confronting the real situation of a finite-energy sum rule involving hadronic data to QCD, 
say a sum rule of the type 

dt p(t) hadronic = F QOD (s 0 ) , (4.59) 

the so values at which the comparison between theory and data should be made are not arbitrary 
p^ |, as is unfortunately often assumed in the literature; optimally, they must be made at the specific 
duality values which are appropriate solutions of the leading OPE constraints relevant to the case one 
is considering. Alternatively, one can choose the highest value of s available by the data (which in 
the real world is fixed by the r mass) so as to minimize the size of the oscillations, and then use the 
MHA to approximate their shape. It would be very interesting to reanalyze r data in these two new 
ways. The case of QCD is of course more complicated than that of our toy model, but our discussion 
of the toy model nevertheless suggests that both ways of analyzing real data should prove useful, and 
possibly lead to an estimate of the errors involved. In addition, one can extend the use of the toy 
model to investigate the MHA for other quantities of interest. Finally, we recall that a slightly more 
sophisticated version of the model does remarkably well phenomenologically 0, indicating that it 
does have something to do with the real solution of large-AC QCD. 
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A APPENDIX 


1. Bernoulli Polynomial 

The Bernoulli polynomials are defined as: 


r —0 


where 


B »(*) =J2( r ) Br a 
T(n +1) 


and i? 2 „ are the Bernoulli numbers : 


r(r + l)T(n-r + l) ’ 
1 


B 0 — 1 , — — , B 2 — — , — — 

0 ’ 2 ’ 6 ’ 30 


1 

42 ’ 


with 


B 2 n+i = 0 , for n > 1. 
Examples of the first few Bernoulli polynomial are: 

B 0 (x) = 1, 

£1(2:) = x - - , 

B 2 {x) = x 2 - x + i , 

6 

£ 3 ( 2 :) = x 3 - ^x 2 + ^x, 

B±(x) = x 4 — 2x 3 + x 2 — . 

ou 


(A.l) 

(A.2) 

(A.3) 

(A.4) 

(A.5) 

(A.6) 

(A.7) 

(A.8) 

(A.9) 


The Bernoulli numbers are related to the Riemann zeta function of even argument: 


B 2n = (-l)" +1 (27r)- 2n 2(2n)!C(2n), for n = 0,1,2,.... (A.10) 

Other important relations for Bernoulli polynomials are: 

B n [x = 0 )=B n , (A.ll) 

B n {x + 1) - B n (x) = nx"” 1 , (A.12) 

Bn{l - X) = (~l) n B n (x ), (A.13) 

B' n {x) = nJ5„_i(x). (A.14) 


It is sometimes useful to recall the trigonometric expansion of Bernoulli polynomials, for N > 1: 


^six) = (-1)™2(2:AT)! ' £ , (A.15) 

V ' n= 1 

iWM = (- 1 )~ +, 2(2JV+l)! Is ri FTT |:^=. (A.16) 

^ 2 n=l 
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